OPERATIONS  RESEARCH 
AND 

INDUSTRIAL  ENGINEERING 


s 


OTIC 

ELECT?? 


AUG  2  5  1980 

A 


COLLEGE  OF  ENGINEERING 
CORNELL  UNIVERSITY 

ITHACA,  NEW  YORK  14853 


fc 


SCHOOL  OF  OPERATIONS  RESEARCH 
AND  INDUSTRIAL  ENGINEERING 
COLLEGE  OF  ENGINEERING 
CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK 


TECHNICAL  REPORT  NO .  440 


January  1080 


INCOMPLETE  BLOCK  DESIGNS  FOR  COMPARING 
TREATMENTS  WITH  A  CONTROL  (IV): 
OPTIMAL  DESIGNS  FOR  p  =  4,  k  =  4 


by 


Robert  E.  Bechhofer 
Cornell  University 

A  jit  C./Tamhane 
Northwestern  University 


/  .. 

Research  supported  by  / 

Army  Research  Office-Durham  Contract' DAAG-: 29- 80-C-0036, 
Office  of  Naval  Research  Contract,'  N0QO14-75-C-0586 
at  Cornell  University  ~ 

and 

NSF  Grant  ENG-77-06112 
at  Northwestern  University 


Approved  for  Public  Release;  Distribution  Unlimited 


TABLE  OF  CONTENTS 


Abstract 

1.  Introduction 

2.  Results  for  p  =  4,  k  =  4 

2.1  List  of  generator  designs 

2.2  Catalog  of  admissible  designs 

2.3  Tables  of  optimal  designs 

3.  Acknowledgment 


References 


ABSTRACT 


The  present  paper  continues  the  study  of  balanced  treatment  incom¬ 
plete  block  ( BTIB)  designs  initiated  in  [1],  [2],  and  [3].  This  class 
of  designs  was  proposed  for  the  problem  of  comparing  simultaneously 
p  >,  2  test  treatments  with  a  control  treatment  when  the  observations 
are  taken  in  blocks  of  common  size  k  <  p  +  1.  A  list  of  generator 
designs,  the  conjectured  minimal  complete  class  of  generator  designs,  a 
catalog  of  admissible  designs,  and  tables  of  optimal  designs  are  given 
for  p  =  4,  k  =  4.  Some  comparisons  are  made  with  admissible  designs 
for  p  =  4,  k  =  3. 


Key  words  and  phrases:  Multiple  comparisons  with  a  control,  balanced 
treatment  incomplete  block  (BTIB)  designs,  admissible  designs,  S-inadmis- 
sible  designs,  C- inadmissible  designs,  minimal  complete  class  of  generator 
designs,  optimal  designs. 
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1.  INTRODUCTION 

The  present  paper  continues  the  study  of  balanced  treatment  incom¬ 
plete  block  ( BTIB )  designs  begun  in  [1],  [2]  and  [3].  This  class  of 
designs  was  proposed  for  the  problem  of  comparing  simultaneously  p  >  2 
test  treatments  with  a  control  treatment  when  the  observations  are 
taken  in  blocks  of  common  size  k  <  p  +  1.  Papers  [l]-[3]  give  the 
background,  motivation  and  notation  for  this  study. 

In  [1]  a  general  theory  of  BTIB  designs  was  developed;  in  [2]  opti¬ 
mal  designs  were  given  for  the  cases  p  =  2,  k  =  2(1)6  and  p  =  3,  k  =  3 
while  in  [3]  optimal  designs  were  given  for  the  cases  p  =  4,  k  =  3  and 
p  =  5,  k  =  3.  In  the  present  paper  we  give  optimal  designs  for  the 
case  p  =  4,  k  =  4;  these  optimal  designs  are  subject  to  the  same  quali¬ 
fication  as  those  given  in  [3] — namely  that  they  are  optimal  relative  to 
the  generator  designs  known  to  us.  However,  we  conjecture  (as  we  did 
for  p  =  3,  k  =  4  and  p=3,k=5  in  [3])  that  we  have  enumerated  all 
of  the  admissible  generator  designs  for  p  =  4,  k  =  4,  and  that  if  ad¬ 
ditional  ones  do  exist  the  incremental  gain  achieved  by  using  the  full 
set  in  place  of  our  set  would  be  very  small. 

In  our  study  of  the  cases  p  =  4,  k  =  3  and  p  =  5,  k  =  3  it  was  neces¬ 
sary  for  us  to  generalize  and  develoo  further  certain  concepts  which  we  had 
introduced  for  the  cases  p  =  2,  k  =  2(1)6  and  p  =  3,  k  =  3.  For  our 
oresent  study  of  the  case  d  =  4,  k  =  4  no  further  generalizations  were  required. 
(See,  however.  Remark  2.3  of  the  present  paper.)  Thus  the  reader  is  referred 
to  [3]  for  the  definitions  of  inadmissibility,  S-inadmissibility  and  C-inadmis- 
sibility  used  in  this  paper.  In  presenting  our  results  for  p  =  4,  k  =  4 
we  hope  to  accomplish  two  objectives:  (a)  To  provide  other  researchers 


in  the  combinatorial  design  area  with  our  list  of  generator  designs  and 
our  conjectured  minimal  complete  class  of  generator  designs  with  the  hope 
that  they  can  supply  additional  generator  designs  (if  any  exist),  and 
more  importantly  that  they  can  propose  a  feasible  method  or  methods  of 
constructing  an  exhaustive  set  of  such  designs,  and  (b)  To  provide  exper¬ 
imenters  with  optimal  (or  nearly  optimal)  designs  that  can  be  implemented 
in  practice. 

The  reader  is  referred  to  Sections  2  and  3  of  [2]  and  Sections  1  and 
2  of  [3]  for  an  exact  statement  of  the  multiple  comparison  problem  under 
consideration,  expressions  for  the  BLUE's  of  the  treatment  effect  differences 
aQ  -  ou  (1  <  i  <  p),  their  variances  and  correlations,  and  an  expression 
for  the  confidence  coefficient  P  associated  with  joint  one-sided  confidence 
interval  estimates  of  the  -  a.  (1  <  i  <  p). 

2.  RESULTS  FOR  p  =  4,  k  =  4 
2.1  List  of  generator  designs 

The  generator  designs  that  we  have  constructed  for  p=4,k=4  (by 
the  methods  described  in  Section  3.2  of  Cl],  or  by  other  methods)  are  lis¬ 
ted  in  Table  2.1.  As  in  [3]  we  have  not  exhibited  equivalent  designs  which 
differ  only  trivially  from  those  given  in  the  tables. 

For  the  generator  designs  in  Table  2.1  we  note  that:  a)  Dg  is 

S- inadmissible  w.r.t.  D,  ,  b)  D.  is  equivalent  to  D„  u  D.  ,  c)  D_  is 

lb  3  4/ 

S- inadmissible  w.r.t.  Dg  and  Dg  u  D^,  d)  Dg  is  S-inadmissible  w.r.t. 

Dg  u  2D4,  e)  Dg  is  S-inadmissible  w.r.t.  u  Dg,  f)  D^  is  S- 
inadmissible  w.r.t.  D^  u  Dg  u  20^,  g)  D^  is  S-inadmissible  w.r.t.  6D4» 
h)  D^,,  is  S-inadmissible  w.r.t.  120^,  and  i)  D  g 


is  S-inadmissible 


Table  2.1 

Generator  Designs  for  p 


Design 


0 


w.r.t.  16D  .  Thus  for  p  =  4,  k  =  4  and  b  ,>  4  it  suffices  to  consider 
unions  of  replications  of  when  seeking  the  optimal  design 

for  a  specified  d/a.  We  conjecture  that  }  is  the  minimal 

complete  class  of  generator  designs  for  p  =  4,  k  =  4. 

2.2  Catalog  of  admissible  designs 

A  catalog  of  admissible  designs  has  been  prepared  based  on  the  set 
of  admissible  generator  designs  given  for  p  =  4,  k  =  4  in  Table  2.1. 

This  catalog  is  given  in  Table  2.2  for  b  =  4(1)27.  It  is  to  be  noted  that 
the  number  of  admissible  designs  increases  rapidly  with  b  for  p  =  4, 
k  =  4;  e.g.,  for  b  =  27  we  have  13  admissible  designs  for  p  =  4,  k  =  4 
whereas  for  b  =  27  we  had  (see  [3])  2  admissible  designs  for  p  =  4,  k  =  3 
and  5  admissible  designs  for  p  =  5,  k  =  3. 

Remark  2.1:  We  note  from  Table  2.2  that  when  appear;,  as  an  admissible 

design  for  a  particular  b,  it  is  always  associated  with  the  smallest  value 
2 

of  t“  and  p  for  that  b;  hence,  the  associated  design  is  always  optimal 
for  that  b  for  d/a  sufficiently  large. 

2.3  Tables  of  optimal  designs 

Optimal  designs  for  p  =  4,  k  =  4  are  given  in  Table  2.3  for 
d/a  =  0.1(0. 1)1.0  and  b  =  4(1)31.  Optimal  designs  that  achieve  a  specified 
confidence  coefficient  1  -  a  are  given  as  a  function  of  d/o  for  1  -  a 
=  0.75,  0.80,  0.85,  0.90,  0.95  and  0.99  in  Table  2.4.  These  designs  were 
found  by  a  complete  computer  search  among  all  admissible  designs. 

Remark  2.2:  We  note  from  Table  2.4  that  for  d/a  -*■  0  and  1  -  a  1 
the  optimal  design  essentially  employs  only  replications  of  D  which 
is  a  BIB  design  among  the  4  test  treatments  augmented  by  a  control  treat- 
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ment  in  each  block.  (An  analogous  phenomenon  was  reported  for  p  =  4, 
k  =  3  and  p  =  5,  k  =  3  in  Remark  3.4  of  [3].) 

Remark  2.3:  As  mentioned  in  Remark  2.1  of  [3]  it  may  be  of  some  interest 
to  compare  designs  with  different  k-values  for  the  same  p-value.  Thus  an 
experimenter  who  is  faced  with  the  choice  of  the  block  size  (subject  to 
the  restriction  that  the  common  block  size  k  <  p  +  1)  may  wish  to  make 
such  broader  comparisons  (which  we  have  not  made  before),  and  rule  out  in¬ 
admissible  designs  using  the  following  generalized  definition:  If 

2  2 

and  D2  are  BTIB  designs  with  parameters  (b  ,k  ,t^ ,p^ )  and  (b2,k2»T2,P2), 

2  2 

respectively,  with  N  =  k  b  <  N2  =  kb,  t1  <  x2»  and  p  £  p? 
with  at  least  one  inequality  strict  then  I)2  is  said  to  be  inadmissible 
w.r.t.  .  (This  definition  is  equivalent  to  the  one  given  in  Remark  2.1 
of  [33  as  a  consequence  of  Theorem  5.1  of  [13.) 

Using  this  definition  it  can  be  verified  that  all  designs  for  p  =  4, 

k  =  3  and  b  =  16  (see  Table  3.3  in  [2])  are  inadmissible  w.r.t.  some  design 

for  p  =  4,  k  =  4  and  b  =  12  (see  Table  2.1  of  the  present  paper),  the 

N  being  equal  to  48  in  both  cases.  Similarly,  all  designs  for  p  =  4,  k  =  3 

and  b  =  20,  24,  28,  32  are  inadmissible  w.r.t.  some  design  for  p  =  4,  k  =  4 

and  b  =  15,  18,  21,  24,  respectively.  We  conjecture  that  the  same  phenomenon 
occurs  for  higher  values  of  b.  This  indicates  that  for  given  p  and 
fixed  kb  =  N,  designs  with  larger  k-values  (k  <  p  +  1)  are  preferable 
since  such  designs  are  "more  complete." 

We  have  obtained  one  example  in  which  a  BTIB  design  for  p  =  4,  k  =  4 
is  inadmissible  (in  the  broader  sense)  w.r.t.  a  BTIB  design  for  p  =  4,  k  =  3. 
(Note:  This  pair  of  designs  was  found  by  comparing  all  admissible  designs 
for  p  =  4,  k  =  4,  b  =  4-75  with  all  admissible  designs  for  p  =  4,  k  =  3, 

b  =  4-130 ,  and  was  the  only  such  design  found.)  If  we  denote  by  D^(3) 


(l<_i<_  5)  and  D^(4)  (1  <*  i  <_  4)  the  designs  for  p  =  4  in  Table  3.1  of 

[3]  and  Table  2.1  of  [4],  respectively,  we  consider  the  designs 


2 

with  b  =  8,  XQ  =  4,  X^  -  4.  Both  designs  have  t  =  4/10,  p  =  1/2 
and  therefore  achieve  the  same  joint  confidence  coefficient.  However,  the 
design  for  k  =  4  has  N  =  32  while  the  design  for  k  =  3  has  N  =  30. 

Thus,  although  each  design  is  admissible  for  its  own  k-value,  the  design 
for  k  =  4  is  inadmissible  w.r.t.  the  design  for  k  =  3.  (This  result  is 
perhaps  not  too  surprising  here  since  the  design  for  k  =  3  is  more  balanced 
than  the  one  for  k  =  4.  in  fact,  the  design  for  k  =  3  is  optimal  for 
d/a  not  too  large  whereas  the  design  for  k  =  4  is  not  optimal  for  any  d/a.) 
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Table  2.2 

Catalog  of  admissible  designs—^  for  p  =  4,  k  =  4 


No. 

of 

blocks 

(b) 

°1 

b^4 

A(1)  -  3 

0  '  3 

=  2 

°2 

b2  *  4 

42)  ■  - 

X<2)  =  0 

" 1 

r» 

"3 

b3  =  6 

A(3)  -  6 

Ao  ■  6 

x<3)  =  i  ! 
1 

°4 

b4  *  1 

44)  - 0 

X<4)  =  1 

Ao 

2 

T 

P 

4 

1 

0 

0 

0 

3 

fl 

0.6061 

0.400 

5 

1 

0 

0 

1 

3 

B 

0.5333 

0.500 

6 

1 

0 

0 

2 

3 

0.4912 

0.571 

0 

0 

1 

0 

6 

n 

0.4667 

0.143 

7 

1 

0 

0 

3 

3 

5 

0.4638 

0.625 

0 

1 

0 

3 

4 

3 

0.4375 

0.429 

0 

0 

1 

1 

6 

2 

0.3810 

0.250 

8 

1 

0 

0 

4 

3 

6 

• 

0.4444 

0.667 

0 

1 

0 

4 

4 

mm 

0.4000 

0.500 

2 

0 

0 

0 

6 

B 

0.3030 

0.400 

9 

1 

0 

0 

5 

3 

H 

0.4301 

0.700 

0 

1 

0 

5 

4 

■ . 

0.3750 

0.556 

2 

0 

0 

1 

6 

D 

0.2821 

0.455 

10 

1 

0 

0 

6 

3 

8 

0.4190 

0.727 

0 

1 

0 

6 

4 

0.3571 

0.600 

2 

0 

0 

2 

6 

0.2667 

0.500 

1 

0 

1 

0 

9  1 

0.2540 

0.250 

11 

1 

0 

0 

7 

3 

0.4103 

0.750 

0 

1 

0 

7 

4 

0.3438 

0.636 

2 

0 

0 

3 

6 

0.2549 

0.538 

1 

1 

0 

3 

mm 

0.2540  | 

0.417 

1 

0 

1 

1 

B 

0.2311 

0.  308 

12 

1 

0 

0 

8 

H 

10 

0.4031 

0.7S3 

0 

1 

0 

8 

m 

8 

0.3333 

O.fcct 

2 

0 

0 

4 

wm 

8 

0.2456 

0.571 

1 

1 

0 

4 

Efl 

6  1 

0.2396 

0.462 

3 

0 

0 

0 

n 

6 

0.2020 

0.400 

—  For  each  number  of  blocks,  the  number  under  D.  (1  <  i  <  4)  in  the  body  of 

the  table  is  the  frequency  f.  with  which  D.  appears  in  the  design 
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Table  2.2  (continued) 


Mo. 

of 

blocks 

(b) 

D1 

bl  =  4 

AoX)  =  3 
=  2 

°2 

b2  =  4 

X<2)  «  4 

0 

X<2)  =  0 

°3 

b3  =  6 

x<3)  .  6 

X<3)  =  1 

°4 

b  =  1 

4 

x<4)  =  0 

0 

x<4>  =  1 

2 

T 

P 

13 

1 

0 

0 

9 

3 

11 

0.3972 

0.786 

0 

1 

0 

9 

4 

9 

0.3250 

0.692 

2 

0 

0 

5 

6 

9 

0.2381 

0.600 

1 

1 

0 

5 

7 

0.2286 

0.500 

3 

0 

0 

1 

7 

0.1922 

0.438 

14 

1 

0 

0 

10 

12 

0.3922 

0.800 

0 

1 

0 

10 

10 

0.3182 

0.714 

2 

0 

0 

6 

10 

0.2319 

0.625 

1 

1 

0 

6 

8 

0.2198 

0.533 

3 

0 

0 

2 

8 

0.1843 

0.471 

2 

0 

1 

0 

12 

5 

0.1771 

0.7OU 

15 

1 

0 

0 

11 

3 

13 

0.3879 

0.813 

0 

1 

0 

11 

4 

11 

0.3125 

0.733 

2 

0 

0 

7 

6 

11 

0.2267 

0.647 

1 

1 

0 

7 

7 

9 

0.2126 

0.563 

3 

0 

0 

3 

9 

9 

0.1778 

0.500 

2 

0 

1 

1 

12 

6 

0.1667 

0.333 

16 

1 

0 

0 

12 

3 

14 

0.3842 

0.824 

0 

1 

0 

12 

4 

12 

0.3077 

0.750 

2 

0 

0 

8 

6 

12 

0.2222 

0.667 

1 

1 

0 

8 

7 

10 

0.2067 

0.568 

3 

0 

0 

4 

9 

10 

0.1723 

0.526 

2 

1 

0 

4 

10 

8 

0.1714 

0.^44 

4 

0 

0 

0 

12 

8 

0.1515 

0.400 

17 

1 

0 

0 

13 

3 

15 

0.3810 

0.833 

0 

1 

0 

13 

4 

13 

0.3036 

0.765 

2 

0 

0 

9 

6 

13 

0.2184 

0.684 

1 

1 

0 

9 

7 

11 

0.2017 

0.611 

3 

0 

0 

5 

9 

11 

0.1677 

0.550 

2 

1 

0 

5 

10 

9 

0.1652 

0.474 

4 

0 

0 

1 

12 

9 

0.1458 

0.429 

Table  2.2  (continued) 


No. 

of 

blocks 

(b) 

D1 

bl  =  4 
X(1)  =  3 

Ao  3 

A^X)  =  2 

°2 

b2  =  4 
X(2)  =  4 

Ao  4 

a<2)  =  0 

°3 

b3  =  6 
A(3)  -  6 

xo  '  6 

X^3)  =  1 

°4 

b4  =  1 
^  =  0 

x<4)  =  1 

*0 

1 

2 

T 

P 

22 

1 

0 

0 

18 

3 

0.3695 

0.870 

0 

1 

0 

18 

4 

0.2895 

0.818 

2 

0 

0 

14 

6 

IS 

0.2051 

0.750 

1 

1 

0 

14 

7 

16 

0.1851 

0.696 

3 

0 

0 

10 

9 

16 

0.1522 

0.640 

2 

1 

0 

10 

10 

14 

0.1455 

0.583 

4 

0 

0 

6 

12 

14 

0.1275 

0.538 

3 

1 

0 

6 

13 

12 

0.1261 

0.480 

5 

0 

0 

2 

15 

12 

0.1143 

0.444 

4 

0 

1 

0 

18 

9 

0.1111 

0.333 

23 

1 

0 

0 

19 

3 

21 

0.3678 

0.875 

0 

1 

0 

19 

4 

19 

0.2875 

0.826 

2 

0 

0 

15 

6 

19 

0.2033 

0.760 

1 

1 

0 

15 

7 

17 

0.1829 

0.708 

3 

0 

0 

11 

9 

17 

0.1501 

0.654 

2 

1 

0 

11 

10 

15 

0.1429 

0.600 

4 

0 

0 

7 

1  12 

15 

0.1250 

0.556 

3 

1 

0 

7 

13 

13 

0.1231 

0.500 

5 

0 

0 

3 

1  15 

13 

0.1114 

0.464 

4 

0 

1 

1 

18 

10 

0.1073 

0.357 

24 

1 

0 

0 

20 

3 

22 

0.3663 

0.380 

0 

1 

0 

20 

4 

20 

0.2857 

0.833 

2 

0 

0 

16 

6 

20 

0.2016 

0.769 

1 

1 

0 

16 

7 

18 

0.1808 

0.720 

3 

0 

0 

12 

9 

18 

0.1481 

0.667 

2 

1 

0 

12 

10 

16 

0.1405 

0.615 

4 

0 

0 

8 

12 

16 

0.1228 

0.571 

3 

1 

0 

8 

13 

14 

0.1204 

0.519 

5 

0 

0 

4 

15 

14 

0.1089 

0.483 

Table  2.2  (continued) 
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D2 

b2  =  4 

X<2)  =  4 

0 

=  0 

°3 

b3  .  6 
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Table  2.4 


1/ 


Optimal  Design— 


to  Achieve  a  Specified  Confidence  Coefficient 
as  a  Function  of  i/a 
for  p  =  4,  k  =  4 


Confidence 

Coefficient 

i/a 

(1-ot) 

0.2 

0.4 

0.6  ! 

0.8 

1.0 

! 

1.2 

1.4 

1.6 

1.8 

2.0 

b=470 

b  =  118 

b=53 

b=30 

b  =  20 

b=14 

b=ll  1 

1 

b  =  8 

b  =  7 

b  =  6 

0.99 

J 

L16 ,0 

28,0 

13,0 

6,0 

5,0 

2,0 

1,0 

2,0 

0,0 

0,0 

1,0 

1,0 

0,1 

1,0 

1  0,0 

1,0 

1,1  i 

0,0 

1  1 

1,0 

b  =  290 

b  =  73 

b=33 

b=19 

b=12 

1 

b=9 

b=7 

b=5 

b=4 

b  =  4 

0.95 

71,0 

19,0 

2,0  1 

3,0 

3,0 

2,0 

0,0 

1,0 

1,0 

1,0 

1,0 

0,1 

0,! 

1,1 

j  0,0 

0,1 

1,1 

0,1  1 

1 

0.0  ! 

0.0 

b=213 

b=54 

3- 

CN 

II 

rQ 

b=14 

1 - 1 

b=9 

b=7 

b=5 

b=4 

b=4 

0.90 

53,0 

13,0 

6,0 

3,0 

2,0 

0,0 

1,0 

1,0 

1,0 

1,0 

0,1 

0,2 

0,0 

0,2 

0,1 

1,1 

0,1 

0,0 

0,0 

0,0 

b  =  168 

b=43 

b  =  20 

b=12 

b=9 

b=5 

b=4 

b  =  4 

b=4 

b=4 

0.85 

42,0 

10,0 

5,0 

3,0 

2,0 

1,0 

1,0 

1.0 

1,0 

1.0 

0,0 

0,3 

0,0 

0,0 

0,0 

0,1 

0,0 

0,0 

0,0 

i  0,0 

b=136 

b  =  35 

b=16 

b=9 

b=6 

b=4 

b  =  4 

b  =  4 

b=4 

b  =  4 

0.80 

34,0 

8,0 

4,0 

2,0 

1,0 

1,0 

1,0 

1,0 

1,0 

1.0 

0,0 

0,3 
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0,2 

0,0 

0,0 

0,0 

0,0 

0.0 

m 

b=28 

b=13 

b=8 

b=5 

b=4 

b  =  4 

la 

b=4 

b=4 

0.75 

n 

7,0 

3,0 

2,0 

1,0 

1,0 

1,0 

1,0 

1,0 

m 

0,1 

0,0 

0,1 

0,0 

0,0 

B 

0,0 

0,0 

— ^The  matrix  in  each  cell  is 


fif2 


where  D  =  u  f.D,  with  b  = 

.,ii 

i  =  l  l 


is  the  optimal  design  for  the  given  value  of  1  -  a  and  i/a. 
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